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' Abstract 

O ' 

We present an electron-identification method using the time-projection chamber and the lead- 
glass calorimeter in the TOPAZ detector system. Using this method we have achieved good electron 
■ identification against hadron backgrounds over a wide momentum range in the hadronic events pro- 

0^ ' duced by both single-photon exchange and two-photon processes. Pion-rejection factors and electron 

efficiencies were 163 and 68.4% for high- Pt electrons and 137 and 42.7% for 1ow-Pt electrons in the 
^ ■ single-photon-exchange process, and 8600 and 36.0% for the two-photon process, respectively. 



1 Introduction 



The TOPAZ detector is located at the TRISTAN e+e" collider of KEK. The center-of-mass energy was 
58 GeV and the integrated luminosity will be more than 300pb _1 by the end of 1994. 

Electron identification at high purity and high acceptance is necessary to tag heavy-quark events. 
In the e + e~ annihilation process, the forward-backward asymmetry of quark-pair production becomes 
maximum at the TRISTAN energy region. A high-accuracy measurement of this was achieved jj], 
In two-photon processes, open-charm production was also studied with high accuracy S. Especially for 
Iow-Pt charm production, electron identification provides the best prove for resolved-photon processes 
, as well as for higher order corrections of the QCD || . 

In order to carry out these studies we need good electron identification over a wide momentum range. 
In particular, the detection of low-momentum electrons is necessary. To meet these requirements we 
have developed a method which uses an energy-loss measurement with a time-projection chamber and 
an energy measurement with a lead-glass calorimeter. In this article we present related details. 



2 The TOPAZ detector system 

The TOPAZ detector is a general-purpose 4-7T spectrometer featuring a time-projection chamber (TPC) 
as its central tracking device H @> ||, f§- A schematic view is shown in Figure |l|. Charged tracks were 
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^internet address: enomoto@kekvax.kek.jp 
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Figure 1: Schematic view of the TOPAZ detector. 
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Components 



Material and size 



Measured Performance 



Magnet 
Beam pipe 
Vertex Chamber 
TCH 
TPC 



TOF 
BDC 

BCL 



ECL 



2.9m</> x 5m, 0.7 X 
10cm(/>, lmmt Be 
Jet-chamber 

Cylindrical Drift Chamber 
2.4m^ x 2.2m 
16 sectors 

10 pad rows, 176 wires / sector 
Ar/CH 4 (90/10) 3.5atm 
Plastic Scintillator 
Streamer tube 

Lead-glass (SF6W), 20 X 
4300 blocks 
| cos 6»| <0.82 
Pb-PWC sandwich 
18 X 

0.85> I cos 01 <0.98 



B=1.0 T 

CT=40/^m 
cr=200Atm 
cr r = 185/im 
cr 2 =335/im 
^£7^=4.6% 



<j Pt /P t = ^(1.5P T ) 2 + 1.6 2 % 

o=220ps 

a=350/im 



a E /E= ^/(8/v^) 2 + 2.5 2 % 
f7 =O.38 o 

a E /E = 6.7% 

(for 26-GeV electron) 

a g =0.7° 



Table 1: Performance of the TOPAZ detector components. 



detected with the TPC placed in an axial magnetic field of 1.0 tesla. The energies of the photons and 
electrons were measured using a barrel lead-glass calorimeter (BCL) and an end-cap Pb-PWC sandwich 
calorimeter (ECL) . A summary of the TOPAZ detector is given in Table @. The coordinate system used 
in this article is as follows: z is defined as the beam-electron direction; x is the vertical axis and y is the 
horizontal axis. The major components of the TOPAZ detector are described in the following sections. 

2.1 The Time-Projection Chamber (TPC) 

The Time-Projection Chamber (TPC) is a three-dimensional tracking chamber which can identify particle 
species by energy-loss measurements. It has a fiducial volume of 2.4m in diameter and 2.2m in axial length, 
comprises 16 multiwire proportional counters (sectors), each equipped with 175 sense wires and 10 rows 
of segmented cathode pads, and is filled with a gas mixture of 90% Ar and 10% CH4 at 3.5 atm fl. 

Charged particles passing through the TPC sensitive region ionize gas molecules and liberate electrons 
along their trajectories. These electrons drift along the direction of the electric field toward the sectors at 
a speed of 5.3 cm//Ltsec, and finally produce avalanches around sense wires; induced signals are detected 
by cathode pads. In order to reduce any space-charge effects, the gating-grid plane is placed above 
the shielding-grid plane. The signals from wires and pads are amplified and are shaped by analogue 
electronics. They are then fed into charge coupled device (CCD) digitizers sampling the wave-forms in 
units of a lOOnsec "bucket" from which the z-positions and energy-loss (dE/dx) information are extracted. 

The spatial and momentum resolutions are studied with cosmic rays as well as Bhabha and 
events. The spatial resolution was obtained to be a xy —185fim and tr z =335/im. The momentum resolution 
was obtained by comparing two measurements in the opposing sectors in cosmic-ray events to be 

<j Pt /P t = v/(1.5P T ) 2 + l-6 2 %. 

The dE/dx resolution was studied using minimum-ionizing pions in the beam events. The obtained 
resolution is a 'dE/dx =4.6% after calibration (discussed later). 
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TOPAZ detector 




Figure 2: Definition of the coordinate system in the TOPAZ detector. 
2.2 The Barrel Lead-Glass Calorimeter (BCL) 

The Barrel Calorimeter (BCL), which comprises 4300 lead-glass Cerenkov counters, has a cylindrical 
structure covering an angular region of | cos#| < 0.82. The calorimeter is divided into 72 modules, 8 in 
the <j> direction and 9 in the z direction, each of which has an array of 60 lead-glass counters 0. 

Each lead-glass counter is made of SF6W and has 20 radiation lengths, so that it can absorb more than 
95% of the shower energies at TRISTAN'S highest energy. The lead-glass blocks are tilted by 1.8° with 
respect to the radial line in the r — <f> plane so that no photons coming from the interaction point escape 
through the counter-to-counter gaps. The Cerenkov light emitted in a lead-glass counter is detected by 
a photomultiplier after passing through a light guide. The signal from the photomultiplier is sent to a 
digitizer and is used for a neutral trigger II. The system is calibrated by using a Xe flash lamp. 

The performance of the BCL has been studied by using Bhabha events. The energy and spatial 

resolution were determined to be a E /E = \J (8/VE) 2 + 2.5 2 % and ere=0.38°, respectively. 

3 Track reconstruction 

The raw data was first corrected for any channel-to-channel variation of electronics. The spatial positions 
of "hits" in the TPC were then determined using pad signals, and particle tracks were searched among 
theses points. Each found track was fitted to a helix, and the charge and momentum were assigned. 

We define the coordinate system in the TOPAZ detector as shown in Figure |[ The global coordinate 
in the TOPAZ detector is defined by x, y and z. The origin of the coordinates is the detector center. 
The £, r] and £ are coordinates fixed on the TPC sector. The £-axis is along the sense-wire direction, the 
f^-axis points in the radial direction at the center line of the sector, and the £-axis is along the beam line. 
The details concerning this section can be found in reference [[HI . 
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3.1 Space-point reconstruction 

The pad signal of each channel was first corrected for any electronics variation, and was then examined 
as a function of the CCD bucket number. Contiguous signals were grouped into a cluster to form a ( hit. 
The peak pulse height and the £ position of the cluster were calculated by fitting the highest three CCD 
buckets to a parabola. Clusters were then searched in the £ direction. If signals in contiguous pads had 
roughly the same ( positions (within 1cm), theses signals were combined into a cluster. 

When a cluster was found, the £ and r\ positions of the cluster were determined. Since the charge 
distribution along the pad row has approximately a Gaussian shape, the £ position was calculated by 
fitting the highest 3 pad pulse heights to a Gaussian. If the cluster had only two pad signals, the pre- 
determined width of the Gaussian was used in the fit. The 77 position of the cluster was calculated by 
taking the centroid of the signals of five wires nearest to the pad row. 

The position resolutions in £ and ( were studied using cosmic-ray tracks; typical values were obtained 
as 0£ = 230/xm and = 340/jm. 

3.2 Track reconstruction 

To find tracks from the reconstructed hits in the TPC, all found space points were first filled in a 2- 
dimensional histogram with <f> and z(R re j) axes such that 

4> = tarT 1 (x I ' \J ' x 1 + y 2 ) 

z(Rref) = ZX(R ref /R), 

where R re f is the radial position of the fixed reference point. The space points associated to a track from 
the origin with a relatively high momentum were clusterized in a bin of the histogram. The track finding 
began by searching for a bin containing more space points than a given threshold. The threshold was set 
at 9 in the beginning, which corresponded to 10 pad raw hits; it was then lowered by one when all of the 
bins above the current threshold were found. Lowering of the threshold was repeated until the threshold 
became two (three pad rows). 

To test whether the space points in the bin formed a track or not, three space points which spanned 
the largest lever arm were chosen and checked as to whether they were on a single helix. If the space 
points were on a helix, other space points in the bin were checked as to whether they were consistent with 
the helix. If not, another set of three space points were chosen and the test was repeated. If a sufficient 
number of space points could be associated they were recognized to form a track and were fitted to a 
helix function. 

Before the helix fitting, the space points were corrected for any dependence on the track angle to the 
sectors. The helix function was parametrized as follows: 

x = —(cos(f)o — cos(cj) + (fio)) + d r cos<f>o + Xq 

K 

y = - (sin<t> - sin{(j> + 4> )) + d r sincj) + Y 

and 

z = —( — )(f)tan\ + d z + Zo, 

where 

(Xq, Yq, Zo) is the position of the pivot, 
k is the inverse of the radius of the helix, 
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</>o is the angle of the pivot to the x axis, 

d r and d z are the distance of the pivot from the true trajectory, and 
A is the the angle which the track makes with the xy plane. 
The parameters to be determined by the fit are 

d r , <j>o, ac, d z , and tanX. 
The fitting of the helix was achieved by minimizing x 2 defined as 

x2 = E(( 6^M )2 + ( £^iM )2) , 

where £j and Zi are the measured position in the i'th pad row and r](r]i) and z(j]i) are the expected 
position by the helix. We checked this \ 2 an d made a final decision as to whether to recognize the space 
points as a track or not. 

After all histogram bins forming tracks were found, the space points which were not recognized as 
tracks were then examined. From the remaining histogram bins, the bin containing the largest number 
of space points was searched for. The inner-most space point in the bin was chosen to be the origin 
(pivot) of the new histogram. The new histogram was defined with wider bins so that the space points 
of a low-momentum track could be clustcrized in a bin as a new track candidate. The bins were tested 
in the same way as described above. 



3.3 Track Refinement 

Although the gating grid prevents the feed-back of positive ions into the drift space and reduces the 
distortion in the electric field in the TPC, a small amount of distortion still remains, resulting in a 
systematic shift in the space-point determination. The dielectric material in the TPC, such as field cages, 
are also a source of distortion. 

In addition to electrostatic distortion, the spatial positions of the TPC sectors may deviate from the 
designed position. This also causes a shift in the determination of the space points. 

The electrostatic distortions in the r</> plane were parametrized as 

= a (r - r ref ) 2 , 

where Ax<p is the position displacement in <j> due to the distortion. The deviations of the sector position 
in the r<p plane (A(, Arj) are parametrized as 

A£ = - Si] 

and 

Ar 1 = d v + 6£. 

The shift in the z direction (dz) was expressed so as to include both the electrostatic distortion and the 
position displacements of the sectors, and was parametrized as 

Az = d z + d v z. 

These parameters were defined for each sector. The parameters of the electrostatic distortions were 
determined based on a comparison of two cosmic-ray data taken when the beam was on and off. The 
parameters of the position displacements of the sectors and the shift in the z-direction were determined 
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based on a comparison of two measurements for a cosmic-ray track by two sectors when the beam was 
on. The space points calculated in the previous sections were corrected for any distortions using these 
parameters. 

Finally, the space points were fitted to a helix again in order to assign the sign of the charge and 
momentum to the track. The momentum in the xy plane (Pr) and along the z-axis (Pl) were expressed 
using the helix parameters as 

Pr = A 

and 

Pl = PrtanX, 

where 

a = — — = 333.56 cmGeV~ l . 
cB 

Q is the charge of the particle, c the velocity of light, and B the strength of the magnetic field (1.0 tesla). 
The sign of the charge was determined from the sign of k. 

The resolution of the momentum measurement was estimated using the cosmic rays jn|. The res- 
olution was studied by comparing two measurements in opposite sectors for a cosmic-ray track. The 
momentum resolution was studied as a function of momentum, and was obtained as 



AP T 



V[(1.5 ± 0.1)P T ] 2 + (1.6 ± 0.3) 2 %. 



The effect of distortion corrections was studied using Bhabha events. The momentum resolution for 
the Bhabha events was measured to be (1.7 ± 0.1)Pt%. 



4 Energy-loss measurement in the TPC 

The energy loss per unit length (dE/dx) is a function of the velocity (/3) and the charge (Q) of a particle. 
In this section we examine the property of the dE/dx curve and discuss what determines the dE/dx 
resolution. In practice, the dE/dx resolution suffers from various systematic shifts due to temperature, 
pressure, high voltage, electronics and so on. We need to improve the resolution as much as possible in 
order to maximize the searchable region. A detailed description of the derivation of dE/dx is given here. 



4.1 Theoretical consideration 
4.1.1 Average energy loss 

When a charged particle traverses a medium it interacts with atoms in the medium. The dominant 
process of energy loss is atomic excitation. The probability of exciting an atom scales as the square of 
the transverse component of the incident particle's electric field with respect to its velocity vector. The 
transition probability, and hence the energy loss, scales as the charge of the incident particle squared. 

Qualitatively, the average energy loss (< dE/dx >) in the material can be described as a function of 
the velocity (/?) and charge (Q) of the particle. The function can be approximately written as 

dE 4™Q 2 e 2 2m e /? 2 7 2 _ 2 

dx m e e 2 (3 2 { n P + (hu>) 2 (3 2 1 2 1 h 

where 

I is the logarithmic mean atomic ionization potential, 
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u) is the plasma frequency of the medium, 

n is the electron density of the medium, 

m e is the electron mass, 

e is the electron charge, and 

Q is the charge of the incident particle. 

The values of < dE/dx > fall as 1//3 2 with increasing the momentum of the particles. The fall off is due to 
the fact that a particle traversing the material spends less time in the electric field of atoms, and, hence, 
transfers less energy as the particle's velocity increases. The ionization becomes minimum at [3j ~ 3, and 
then < dE/dx > rises logarithmically as the momentum increases. This "relativistic rise" is the result 
of the relativistic increase in the transverse electric field of the incident particle. The < dE/dx > finally 
reached the plateau due to the "density effect" caused by the polarization of the medium. 

4.1.2 Energy-loss distribution 

There are two contributions to the shape of the energy-loss distribution: atomic excitations and 8 — 
ray production. The 5 — ray contribution falls inversely proportional to the energy-transfer squared 
(l/(AE) 2 ), because this is the classical Rutherford scattering for free electrons. This process makes a 
long tail, called Landau tail, in the energy- loss distribution up to the kinematical limit. The peak of 
energy loss distribution is due to atomic excitation, which is the dominant energy-loss process. Argon, 
which is the main component of the TPC gas, has three shell structures: K, L, and M. The energy 
levels are 3.20 keV, 248 eV, and 16 - 52 eV, respectively. Since the collision probability with an electron 
is roughly proportional to its energy level, the lower energy-level excitations give a narrower dE/dx 
distribution. 

4.2 dE/dx measurement 

4.2.1 Derivation of the energy loss 

The energy loss was measured using wire signals of the TPC. First, the wire signal clusters were searched 
in the z direction in the same way as the pad clusters. The pulse height of a cluster was determined 
from the peak amplitude of a parabola fitted to the highest 3 CCD buckets in the cluster. The z position 
of the cluster was also determined from the fit. The cluster was then associated to a track using the z 
information. If the cluster was within 1 cm of the track in the z direction, and not associated to other 
tracks, the cluster was linked. 

The pulse-height amplitude was converted into the energy-loss value by two steps. The first was a 
correction for any non-linearity of the electronics. The second was a correction for the wire-to-wire gain 
variation. Details concerning the calibration are described in the following subsection. 

Since a TPC sector has 175 sense wires, the energy loss for a track was sampled every 4-mm segment 
up to 175 independent data. A single variable, which is similar to the most probable energy loss, was 
deduced from the data sample. The final dE/dx was then derived by making various corrections, which 
are discussed later. 

4.2.2 Calibrations 

Prior to the experiment, after the electronics were calibrated by pulsing the shielding grid wires with 
various test-pulse amplitudes, the calibration curve for each channel was determined. The linearity of 
our electronics is better than 0.5% below the saturation point. The remaining non-linearity, appreciable 
near to the saturation point, is corrected using this calibration curve. 
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Figure 3: Pulse-height distribution obtained by the TPC wire signals for minimum ionizing pions. 

After an electronics calibration, a correction for the wire-to- wire gain variation is made. Each sector is 
equipped with three rods having a 55 Fe X-ray source for each wire. Each rod can be moved pneumatically 
behind the hole on the cathode plane to irradiate wires. The wire-to-wire gain variation is corrected by 
calibrating the pulse height for the main peak of 5.9 keV for each wire. This calibration reduces the gain 
variation over a sector within a 3% level, and also gives the factor necessary to obtain the absolute value 
of the measured dE/dx. 

4.2.3 The 65%-truncated mean 

Figure |3| shows the measured energy-loss distribution in the TPC for minimum ionizing pions. The shape 
is asymmetric and has the long tail at the higher side. Therefore, a simple mean of the dE/dx sample for 
a track is not a good parameter for particle identification. The most probable energy loss or the mean 
value around the peak is a better parameter. We adapted the mean of the lowest 65% in the data sample 
(the 65%-truncated mean) as an estimator for the energy loss of a track. In the following we use dE/dx 
as in this meaning. 

4.2.4 Corrections 

To achieve a good resolution in dE/dx measurements, one must reduce the systematic shifts of the gain. 
In this section we discuss the systematic shifts of dE/dx and how we correct them. 

In the TPC, measuring the energy loss is an indirect process. An implicit assumption is that the 
number of ionization electrons produced per track length is proportional to dE/dx of the charged particle. 
As electrons in the gas drift to the sectors they defuse; some fraction is absorbed by impurities, such 
as (?2j in the gas. Since the maximum drift length is more than 1 m, this effect is not negligible. We 
determined the attenuation factor using minimum ionizing pions in the momentum range between 0.5 
and 0.6 GeV/c. Figure |J shows the measured value of dE/dx as a function of the drift length. We correct 
this effect using following function: 

P H correct — PUmeasured X (1 4~ Cli/), 

where 
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Figure 4: dE/dx as a function of the drift length. 



PH is the pulse height, and 

L is the drift length. 

The attenuation factor (Ci) is typically 7% per meter. 

The gain shift due to the sample thickness dependence was also corrected using the minimum ionizing 
pions. Figure |E| shows dE/dx vs. the path length on the logarithm scale. This systematics was corrected 
using 

P H correct — PPmeasur cd x (1 - C2 x log(X)), 

where 

X is the path length normalized by 4 mm. 

The factor Ci was determined to be 0.2. 

The electric fields near to the sense wires also affect the avalanche process. The dependence of the 
pulse height on the high voltage was parametrized as 

SV 

P P correct — P Hmeasured X (1 C3 X ) , 

where 

Vo is the nominal high voltage (1970V) on the sense wire, and 

SV is the variation of the voltage. 

The correction factor (C3) of the pulse height was determined using data taken in low high- voltage 
cosmic-ray runs. A 1% high-voltage change corresponds to a 17% change in the gain. 

The electrons approaching the sense wires pass through the avalanche process. As the gas density 
increases the mean-free path of electrons becomes shorter and the soft collisions in which the electrons 
do not ionize the atoms increase; hence, the amplification factor drops. The gas density is proportional 
to the pressure and inversely proportional to the temperature. The pressure was well controlled, and its 
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Figure 5: dE/dx as a function of the sample thickness in logarithm after a segment-length correction. 



5 
4.8 
4.6 
4.4 
4.2 

4 
3.8 
3.6 
3.4 
3.2 

3 




-0.8 



-0.6 -0.4 -0.2 0.2 
Deviation of density (%) 
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variation was less than 0.04% throughout this experiment. The temperature of the gas was increased by 
about 2°C during the summer. Figure ^| shows dE/dx vs. the variation in the gas density. From the 
temperature change, we found that 



where 

po is the nominal gas density, and 
Sp is the variation of the gas density. 

The C4 value was 6.0. A 1°C temperature change corresponds to about a 2% change in the gain. 

After these corrections, we again carried out a wire-to-wire gain calibration. The event sample with 
minimum ionization pions were used for this purpose. After the pulse-height distributions for each 
wire were made and the correction factors were derived, a run-dependent correction was carried out. 
There were still unknown factors which changed the energy-loss values. Each step used a single run 
corresponding to a few 1000 beam events. 

Finally, we derived the dE/dx distribution for momentum between 0.5 and 0.6 GeV, which was the 
minimum ionizing pion region. There were two peaks corresponding to pions and electrons. For several 
reasons, the relationship between these two peak values changed nonlinearly. One reason was an over- 
value setting of the threshold values in the electronics during some experiments. The other was a selection 
bias in the event pre-selection code in order to reduce the data size, mainly due to the vacuum condition. 
These caused a nonlinear offset in the pion dE/dx value. We thus corrected these by normalizing them 
to the initial experimental values in which the threshold values were considered to be sufficiently low. 
Furthermore, we fitted the width of two dE/dx peaks using the common ratio to the center values of the 
peaks using two Gaussians. The fitted width was considered to be the effect of this nonlinearity; the % 2 s 
calculation is described later. 

4.2.5 Velocity dependence of dE/dx 

The dE/dx curve as a function of the particle velocity was determined from the data. Although theo- 
retical calculations of the dE/dx curve reproduce the data fairly well, they are not adequate for particle- 
identification purposes, especially in the relativistic rise region. The data used for the determination were 
those from protons for the small-/3 region, cosmic-ray muons and minimum ionizing pions for low f3j, and 
electrons for the relativistic rise region. The data are fitted to a ninth-order polynomial parametrized as 



where the a are given in Table Figure [| shows the dE/dx vs. momentum for tracks in more-than- 
three-track events. The bands of electrons, pions, kaons and protons can be clearly seen. 

4.2.6 dE/dx resolution 

The resolution of dE/dx was studied using the minimum ionizing pions. Figure |§| shows the dE/dx 
distribution of pions with momenta between 0.5 and 0.6 GeV/c. The resolution was estimated by fitting 
a Gaussian function to the distribution, and was obtained as u dE / dx = 4.6%. The resolution for Bhabha 
events was also measured as a dE/dx = 

4.4%. 



PH, 



correct — 



PH measure( i X (1 + C4 X ), 
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coefficient value 



co 


6.2133 


Cl 


2.1340 


C2 


-1.6116 


C3 


-0.19783 


c 4 


0.0363802 


C5 


-0.35785 


c 6 


1.4957 


C7 


-0.16234 


C8 


-0.95077 


eg 


0.379 



Table 2: Coeficients of the 9'th-order polynomial function. 
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Figure 7: dE/dx vs. momentum for tracks in more-than-one-track events. 
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Figure 8: dE/dx distribution for minimum ionizing pions. 



5 Energy measurement in the BCL 
5.1 Clustering 

The electromagnetic clusters detected with the BCL were used for electron identification. The method 
for electromagnetic energy clustering is described as follows: 

1. We searched the maximum energy counter among 8 neighboring counters (at the +z or -z edge of 
the BCL, 5 counters) for counter by counter; 

2. If the neighboring maximum counter energy was higher than its energy, the counter was merged to 
the neighboring one; 

3. After the above procedures were iterated, the clusters were formed; 

4. The cluster energy was calculated by summing the total counter energy in the cluster; 

5. The azimuth angle ((f)) and the z-position were calculated using the energy- weighted mean of the 
counter positions; 

6. The depth of the cluster center was calculated according to X = ln(E/E c ) + 1.5 (radiation length), 
where E is the cluster energy and E c (12.6 MeV) is the critical energy of the BCL (lead glass); 

7. The radial position (R) of the cluster was calculated according to the formula R — Rbcl + X ■ 
Xcos9, where Rbcl and Xq are the inner radius and the unit radiation length of the BCL, i.e. 
176cm and 1.7cm, respectively. Here, 9 was assumed to be that of the maximum energy counter in 
the cluster, at first; 

8. From z and this R, the was re-calculated; and 

9. R and 9 calculations were iterated two times. 
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5.2 Energy calibration 

The monitoring system for the gains of 4300 counters comprises just one light source and an optical-fiber 
system. As the light source, a xenon flash lamp having a 3mm gap was chosen, because the emission 
spectrum resembles Cerenkov light. The amount of light received at a counter corresponds to about 7 
GeV in energy. In order to correct for any variation in the light intensity of the xenon flash lamp, reference 
counters were placed in a temperature-controlled box. At each reference counter a NaI(Tl)- 241 Am light 
pulser was attached to measure the gain independently. The gain variation can be corrected to the 0.5% 
level. 

The gain of each counter was calibrated with 4, 2, and 1 GeV electron beams from the IT1 beam line 
jL§ at the 8-GeV Accumulator Ring (AR) at KEK. 

We have accumulated Bhabha scattering events, after standard cuts of requiring two clusters of 
deposited energy greater than 1/3 of the beam energy at | cos0| < 0.77 satisfying an acolinearity angle 
being less than 10°. The obtained resolution was 4.5%, apparently dominated by the nonuniformity in 
the gains of individual counters. The angular resolution was obtained to be 0.36°. Calibration of the 
individual counters was carried out using these Bhabha events. However, not all of the counters have 
been calibrated yet, due to a lack of statistics. 

The electron energy-loss in the material in front of the BCL was investigated by studying the E/P 
for electrons. In this study we used an experimental two-track sample selected by the following criteria, 
expecting these to be dominantly e + e~ — > e + e~e + e~ events: 

1. The number of charged tracks with Pt > 0.15 GeV, the polar angle | cos6*| < 0.83, and the number 
of degrees of freedom (N.D.F.) of charged-track fitting > 3 had to be 2; 

2. The number of the BCL clusters with E c i uster > 100 MeV had to be at least 1; and 

3. The visible energy (E V i S ) of the event had to satisfy E V i S < 30 GeV, where both the charged tracks 
in the TPC and the clusters in the BCL were used. 

Most of the selected events were expected to be two-photon events of (e + e~)e + e - , (e + e~)/i + /i~, and 
(e + e _ )7r + 7r~. The remaining ones were u + u~ , t + t~, e + e~j, etc. The used data corresponded to an 
integrated luminosity of 24.5 pb" 1 at y/s = 50 ~ 60.8 GeV. These selection cuts left 9120 two-track 
events. 

To separate electron tracks from muons and pions, we used energy- loss (dE/dx) information from the 
TPC and TOF as follows: (1) The electron-candidate tracks were required to have dE/dx in the range 
of 5.5 < dE/dx < 7.5 keV/cm; (2) the confidence level (CL) for the electron hypothesis by the TOF 
information had to be > 0.01. We then calculated the E/P ratio of the energy measured by the BCL 
and the momentum measured by the TPC. The TPC-BCL combination was chosen when the distance 
between the BCL cluster and the extrapolated track was shortest. We only calculated the E/P for the 
TPC-BCL combination whose distance was less than 10cm. 

To determine the BCL energy correction we fitted the E/P distribution with a second-order polynomial 
of the BCL energy. We used electron candidates which were in the energy range E < 4.0 GeV and the 
polar angle range |cos#| <0.77. The E/P-versus-energy distribution with the fitting curve is shown in 
Figure ||. Since real electrons should concentrate on E/P = 1, we obtained the following expression from 
the fitting result: 

^ = 0.430 - 0.354 • E + 0.072 • E 2 for E < 2.3GeV, and 
= for E > 2.3GeV. 

The E/P distributions in the momentum range 0.4 < P < 4.0 GeV, E/P and momentum scatter plots 
with and without the BCL correction are shown in Figure [Io| 
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Figure 9: E/P ratio as a function of the BCL energy (filled circle). The solid line was obtained from a 
polynomial fitting (described in the text). 



6 Detector simulation 
6.1 Event generation 

For single-photon-exchange process, we used JETSET6.3 as an event generator p3| . The parameters 
used in this program were fitted for general event shapes |l4| . For two-photon processes, we used own 
event generator of which the details were described in reference p[ [t5| [i6[. 



6.2 Monte-Carlo simulation 

The particles generated by the above-mentioned event-generation program were traced through the 
TOPAZ detector simulator. This detector simulation was designed to take into account all conceivable 
interactions which the particles might experience in the detector, and to mimic the detector response 
to them as closely as possible. The structures and materials of the detectors are precisely coded in the 
program and each particle is propagated by a small step from the interaction point. At each step, such 
interactions as decays, multiple scatterings, and bremsstrahlungs take place according to the probabilities 
associated with them. 

In the TPC region the averaged energy loss in the gas is calculated according to a formula determined 
from the minimum ionizing pions in real data, knowing the velocity and charge of the particle. The 
Landau fluctuation is taken into account by smearing the energy loss so as to be consistent with the 
experimental distribution. The smeared energy loss is translated to a pulse height on the corresponding 
wire to be sampled in order to mimic the digitized outputs. Because the time profile of the wire signal 
was experimentally found to be independent of the track angle, a single pulse-height shape function is 
used to calculate the pulse height in each time bucket. The pulse height is digitized if it is above the 
threshold, and is recorded as wire data. Pad signals are generated using the pulse heights of the relevant 
five wires near to the pad rows. The avalanches on the five wires induced a charge along the pad row. 
We assume that the induced charge is a Gaussian distribution, and is proportional to the corresponding 
wire pulse heights. The width of the Gaussian and the weight factors for the five wires are determined 
experimentally fTij] . 
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Figure 10: E/P distributions in the momentum range 0.4 < P < 4.0 GcV [(a)-l and 2] and Scatter 
plots of E/P vs. energy [(b)-l and 2]. Indices 1 and 2 denote that without and with the BCL correction, 
respectively. The equation of the correction function is given in the text. 
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To simulate the electromagnetic shower development in the material, the EGS code is widely used 
JL8] ] . The cut-off energy in the EGS program, when it is used, must be set at a value sufficiently lower 
than a critical energy for an accurate simulation, which consequently needs more CPU time. In the 
BCL, the electromagnetic shower is simulated by using the Bootstrap method to save CPU time fll9fl . 
This method utilizes pre-generated showers to obtain the shower density in cells of suitable size in the 
calorimeter. The development of a shower with an energy lower than the cut-off is replaced by them. 
Thus, once many sets of the frozen showers are prepared for various energies, the shower simulation can 
be carried out as accurately as that of the EGS code, while at the same time considerably reducing the 
CPU time. 

For a nuclear interaction GEISHA 7 is used p0|| . Short-lived particles such as K s , ir°, and A are 
decayed in the simulator program while keeping secondary vertex information. 



7 Hadronic-event sample 
7.1 Trigger 

The charged trigger required at least two charged tracks with Pt > 0.3-0.7 GeV and an opening angle > 
45-70°, depending on the beam conditions. The neutral energy trigger required that the energy deposit 
in the BCL had to be greater than 2-4 GeV, or that in the ECL had to be greater than 10 GeV. Details 
concerning the trigger system can be found in reference 0. 



7.2 Single-photon-exchange process 

Our selection criteria for the single-photon-exchange process are as follows: (1) the number of good tracks 
had to be > 5, where a good track was required to have (i) a transverse momentum to the beam axis 
larger than 150 MeV, (ii) a distance of closest approach to the beam axis less than 5cm in the xy-plane, 
as well as in the z-direction, and (iii) a polar angle satisfying 0.02 < | cos6>| < 0.77 to be well contained 
in the detector; (2) the visible energy (E V i S ) had to be > the beam energy (E oeam ), where E V i S is the 
sum of the momenta of the good tracks and the energies of BCL and ECL clusters of energy greater than 
100 MeV; (3) the longitudinal momentum balance had to be | EP 2 1 /Ebeam < 0.4, where the P^s were 
longitudinal momenta of the good tracks or longitudinal energies of the BCL and ECL clusters; (4) the 
larger of the hemisphere-invariant masses had to be greater than 2.5 GeV, where the two hemispheres 
were divided by a plane perpendicular to the thrust axis; and (5) those events with two or more large 
energy clusters of greater than 0.5Eb eam were discarded. In total, 23,783 hadronic events were obtained. 
The corresponding luminosity was 197 pb _1 and the average center-of-mass energy was 58 GeV. The 
details can be found in reference [|[ pl| . 



7.3 Two-photon process 

The hadronic events produced by two-photon processes were selected based on the following criteria: (1) 
The number of good tracks had to be at least 4, where | cos(9| should be less than 0.83; (2) the position 
of the origin of the event (i.e., the event vertex), reconstructed from all tracks, had to be within 1.5cm 
in the xy-plane and within ±2. 0cm in z-direction; (3) E V i S had to be less than 30GeV, where only the 
BCL clusters were used for neutral energy; (4) the mass of the hadronic system (W V i S ) had to be greater 
than 3GeV; and (5) the energy of the most energetic cluster appearing in the BCL, ECL, and FCL had 
to be less than 0AEb eam (anti-tag condition). These selection cuts left 23,779 events at an integrated 
luminosity of 95.3pb _1 . The average center-of-mass energy was 58 GeV. The details can be found in 
reference [pi. 
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Figure 11: Distributions of the closest approach in the Monte-Carlo simulation; the solid histogram is 
for prompt electrons and the dashed one is for pair conversion and Dalitz decays. 

8 Electron selection for a single-photon exchange 
8.1 Good track selection 

The electron tracks which satisfied good-track criteria were selected by the dE/dx information from the 
TPC and the E/P ratio measured with the TPC and BCL in a broad momentum region. 
The selection criteria for good tracks were followings: 

1. The closest approach to the beam axis in the xy-plane must be less than 1.0cm; 

2. The closest approach to the interaction point in the z-direction must be less than 4.0cm; 

3. The absolute value of the cos 9 of the track must be within between 0.02 and 0.83, where 9 is the 
polar angle of the tracks with respect to the initial electron direction; 

4. The transverse momentum with respect to the beam axis must be greater than 0.15 GeV/c; 

5. The number of degrees of freedom (N.D.F.) of the track fitting must be greater than 3, with a 
maximum of 15; 

6. The number of wires for the dE/dx calculation (65-% truncated mean) must be greater than 30, 
with a maximum of 114; and 

7. The momentum must be greater than 0.8 GeV/c. 

Conditions (1),(2), and (5) require good measurements of the momentum of the tracks. (1) also 
reduces the background. Condition (3) requires the tracks to be within the TPC active volume. (4) 
requires that the tracks are not curled up in the TPC volume. (6) is the condition for the sufficient 
resolution of dE/dx. (7) is added to reduce the background. 

Figure [n] shows the difference of the closest approach to the beam axis in the xy-plane for prompt 
electrons and the background tracks in a Monte-Carlo simulation. Figure [l2| shows the difference in the 
momentum distribution for prompt electrons and the background tracks in a Monte-Carlo simulation. 
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Figure 12: Distributions of the momentum in the Monte-Carlo simulation; the solid histogram is for 
prompt electrons and the dashed one is for the background. 

After the tracks were subsequently extrapolated to the radius of each BCL cluster center, we calculated 
the distance to the clusters. The cluster was matched with the track having the closest distance. E/P 
was calculated using a matched track and cluster. 



8.2 Electron selection 

The electron selection was carried out based on two criteria. 

First, we calculated the %|s for each hypothetical particle using a dE j dx information and a momentum 
measurement from the TPC, where i indicates an electron, muon,pion, kaon, or proton. The definition 
of Xi is 
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dx meas. " 171 

dE 



meas.: respectively. They 

include inefficiency measurements in the TPC. A^ also includes the effect of the nonlinear shift 

J dx meas. 

described before. The most probably x\ is given by minimizing it while modifying f3j. The expected 
dE/dx is given by a polynomial function of (3-f. The details are described in Section 4.2.5. 

The electron-candidate tracks were selected by requiring that xl must be less than 3.0, where the 
number of degrees of freedom (N.D.F.) is 1. Figure |l3| shows the E/P distributions both before and after 
this selection. The peak of the electrons is enhanced. 

The rejected tracks in this selection (di^/dx-hadron) are used for hadron background estimations. 

Second, we required the electron-candidate tracks to have a distance of less than 5.0cm with respect 
to the calorimeter cluster positions, and calculated the cluster width, which is the energy-weighted r.m.s. 
of counter positions with respect to the extrapolated track position. We also carried out selection for 
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Figure 13: E/P distribution of the experimental data; the dashed histogram is for the before dE/dx 
selection and the solid one is for the after dE/dx selection. 



the cluster width, i.e., we selected the clusters having widths of between 1.0cm and 10.0cm. These 
distributions for electron-candidate tracks and hadron tracks are shown in Figure 
The E /P distributions of before and after these selections are shown in Figure 

8.3 Background estimation 

8.3.1 7-conversion rejection 

Since the largest background source was electrons from ~ ; -conversions at the material in front of the TPC. 
we rejected the dominant part of such electrons by the following methods. 

We reconstructed secondary vertices (V 0, s) with an electron-candidate track and one of the unlike sign 
charged tracks, and then calculated the invariant mass of each V°, assuming that its daughter particles 
are electrons. For the ^-reconstruction, two kinds of vertices, i.e., non-crossing and crossing cases in 
the xy-plane (perpendicular to the beam axis), were searched. In the former case, the combinations 
with a distance in the xy-plane of less than 4.0cm and a z-distance of less than 2.0cm were selected; we 
required that the deflection angle between the momentum vector of two tracks and the position vector of 
the closest point be less than 5.0 degrees. In the latter case, we selected one vertex which had a smaller 
deflection angle, and required that the z-distance be less than 2.0cm and the deflection angle be less than 
5.0 degrees. 

We then rejected the tracks in the pair if its invariant mass was < 50 MeV in the former or < 200 
MeV in the latter cases, respectively. 

These distributions are shown in Figure [l6|. 

8.3.2 Hadron background estimation 

The hadron background on the E/P plot was estimated using the dE / dx-hadions. The distributions 
were normalized using the entries between 0.0 and 0.64. The normalized E/P distributions of the dE/dx- 
hadrons were subtracted from the E/P distributions of the selected electron candidates. The numbers of 
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Figure 14: Distributions of (a)the minimum ce between the track-extrapolated point an 

cluster, and (b)the shower width with respect to the track; the solid histograms and the dashed histograms 
are the total tracks and the electron-enhanced samples of the hadronic events, respectively. 
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Figure 15: E/P distribution of the experimental data; the solid and the dashed histograms are after and 
before the shower selection, respectively. The dotted one is for the rejected tracks by this selection. 
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Figure 16: Distributions of the secondary vertices in the Monte-Carlo simulation; (a)the distance in the 
xy-plane of the not intersected case, (b)the distance in the z-direction of the not intersected case, (c)the 
distance in the z-direction of the intersected case, (d)the deflection angle of both cases, (e)the invariant 
mass of the not-intersected case, and (f)the invariant mass of the intersected case. 
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Figure 17: Example of the hadron background estimation; the solid and dashed histograms are the 
selected electron candidates and the dE/dx-hadrons which is described in the text, respectively. 

electrons were obtained by counting the entries between 0.74 and 2.0 in the subtracted E/P distributions. 
An example of the hadron-background estimation is shown in Figure [I?]. 



8.3.3 Remaining 7-conversion estimation 

There remained electrons resulting from 7-conversions which escaped from V° reconstruction. They were 
presumably energy- unbalanced pairs, for each of which the lower Pt track was not reconstructed by the 
TPC. We estimated the failure rate of the V° reconstruction by a Monte-Carlo simulation. The failure 
rate, 77 = N^ f ailure /N^ f econstructed , was estimated for each P T -bin. In the calculation, we included 
contributions from Dalitz decay(7r° — ► e + e~j) as well as the conversion electron pairs. Since we did not 
apply any vertex position cut in the conversion-pair reconstruction, this factor included the effect on the 
Dalitz pairs. In order to estimate the number of remaining background tracks from 7-conversions, we 
multiplied the number of reconstructed conversion pairs in the experiment by r\ in each Pt bin. 



8.4 Pion-rejection factor, electron efficiency 

The electron selection in the single-photon-exchange process is summarized here. The pion-rejection 
factors were derived using single-track multi-hadron events generated by Monte-Carlo simulations. In 
order to derive the efficiency and the factor, we required E/P to be within 0.72 to 2.0 for the selected 
tracks. All of the tracks which passed the above criteria were counted. The definition of the selection 
efficiency is the ratio of the number of electrons counted (as above) to the number of charged tracks 
before all of the selection criteria. The pion-rejection factor is defined as the ratio of the number of pions 
before selection to the number of the remaining pions with an E/P of between 0.72 and 2.00. 

This value for single-track events is shown in Figure together with the electron efficiency. Since 
the dE/dx curve of the pions is close to that of the electrons at the higher momentum region, the pions 
satisfy the dE/dx selection criterion. Thus, the pion rejection factor is worse in the higher momentum 
region. The average values were 80 for pion rejection and 68.4% for electron efficiency. In order extract 
these two averages, we assumed the momentum spectrum in the multi-hadron events. 
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Figure 18: Pion-rejection factor and electron efficiency for single-track events; the solid curve is the 
pion-rejection factor (left scale) and the dashed one is the electron efficiency (right scale). 



In practice, further rejection was carried out. The pions remaining after the E/P selection were 
estimated using the tracks rejected by the dE/dx selection (AE/cfe-hadrons). The E/P distributions 
of the AE/cix-hadrons were normalized using the entries in the side-band (0.0 <E/P< 0.64); then the 
distribution from the E/P distribution of the selected tracks was subtracted. This procedure improves 
the pion rejection factor by more than factor 3. 

For multi-hadron events, we show these values for two cases, i.e., low-Py (Pp < 0.8 GeV) and high-Py 
(Pr > 0.8 GeV). They are shown as a function of the electron momentum in Figures 19 and|2fj|. The same 
E/P selection was carried out as that described above to derive these Figures. For the higher momentum 
tracks, the pion-rejection factors are rather worse because of the dE/dx effect (as described before). To 
summarize, the pion-rejection factors of high-Pp tracks and low-Pr tracks were estimated to be 164 and 
137 on the average, where the electron efficiencies were 68.4 and 42.7%, respectively. Practically, we 
also used additional pion rejection by fitting (as mentioned above). We can thus we can expect better 
rejection of the pions. 



9 Electron identification in two-photon processes 
9.1 Electron selection 

In two-photon processes electron tracks were identified by combining the information from the E/P ratio, 
dE/dx, and the TOF as follows: 

1. Charged tracks in the TPC were extrapolated to the BCL. We then selected for each TPC track 
the BCL cluster which was the closest. The E/P of each of the so-selected TPC-BCL combinations 
had to satisfy 0.75 < E/P < 1.25. 

2. dE/dx was calculated for the electron-track candidate to be used in the subsequent electron count- 
ing. The resultant dE/dx was mostly contained in the range 5.5 < dE/dx < 7.5 keV/cm. 

3. We used information from the TOF and required the track to have a confidence level(CL) of 0.01 
or better for the electron hypothesis. 
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Figure 19: Pion-rejection factor and electron efficiency for \ow-Pt tracks; the solid curve is the pion- 
rejection factor (left scale) and the dashed one is the electron efficiency (right scale). 
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Figure 20: Pion-rejection factor and electron efficiency for high-Pr tracks; the solid curve is the pion- 
rejection factor (left scale) and the dashed one is the electron efficiency (right scale). 
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The energy-loss (dE/dx) information from the TPC enables us to separate electrons from hadrons in 
the momentum region(Pr < 3GeV). The E/P ratio of the energy (E) measured by the BCL and the 
momentum(P) measured by the TPC can separate electrons clearly above 0.4 GeV. The TOF was useful 
to resolve electrons from kaons and protons in the overlapping region of dE/dx. 

The performance of electron identification in the momentum range 0.4 < Pt < 3.0 GeV and the 
polar-angle range | cos 6 1 1 < 0.77 is demonstrated in Figures |l], where various distributions (the closest 
TPC-BCL distance(Figure |l|-(a)), thc E / p ratiojTigure |l]-(b)), the CL for the electron hypothesis in 
the TOF(Figure |l|-(c)), and the dE/dx (Figure |lf- (d))) are shown. In these figures, electrons coming 
from 7-conversions were already rejected. The methods for 7-conversion rejection are described in the 
next section. Notice that these figures were obtained with the electron candidates selected by all of the 
cuts except for that on each plotted quantity. Notice in particular that Figures ^l)-(a) to (c) were made 
with an dE/dx cut (5.5 < dE/dx < 7.5 keV/cm) which was used only for the purpose of displaying. 
As can be seen from the dE/dx distribution in Figure pl"r(d), two peaks corresponding to electrons and 
pions are clearly separated. 

We counted the numbers of electrons in each Pt bin by fitting the dE/dx distributions with double 
Gaussians bin by bin. The Pt binning was selected so as to approximately equalize the number of entries 
in each bin. The dE/dx distributions are shown in Figures ^2|. 



9.2 Reject 7-conversions 

The electrons from the 7-conversions were rejected by secondary vertex (V°) reconstruction. The V°'s 
were reconstructed from all combinations of two tracks, and the invariant mass of each V° were calculated. 
When two tracks did not intersect at the closest point in the xy-plane, the distances of the two tracks 
at the minimum-distance position were required to be less than 7 cm in the xy-plane and 3 cm in the 
z-direction. When two tracks intersected, we chose from the two crossing points that with the shorter 
z-difference, and required it to be less than 1.5 cm. We then rejected the tracks in the pair if its invariant 
mass was < 80 MeV in the former or < 150 MeV in the latter cases. 

We also required the closest approach of each electron-track candidate to the event vertex in the 
xy-plane to be < 0.5-1.5 cm, depending on Pt- 



9.3 Remaining 7-conversion estimation 



The remaining 7 conversions and Dalitz decays were estimated using the failure rate of the V° recon- 
struction which was obtained by a Monte-Carlo simulation. The failure rate of the V° reconstructions 
(described in section 8.3.3 ) was estimated for each Pr-bin, and was typically ^0.4. The number of remain- 



ing background tracks was estimated by this ratio multiplied by number of the reconstructed conversion 
pairs in the experiment in each Pt bin. The 7-conversion and Dalitz decay background mainly occupy 
the low-Pr region. 



9.4 Pion rejection factor and electron efficiency 

The electron selection in the two-photon processes is summarized here. The electron efficiencies for 
single-track events and for multi-hadron events which were generated by the Monte-Carlo simulation are 
shown in Figures ^3] and ^4|, respectively, as a function of the transverse momentum of the track. The 
average values of the efficiencies were 34.2% for single-track events and 36.0% for multi-hadron events. 
The former value was obtained assuming the Pt distribution in the multi-hadron events. 

The pion-rejection factors were also derived with single-track and multi-hadron events. For single- 
track events we obtained a pion-rejection factor of 15000 for the high momentum region of 1.5 < Pt < 3.0 
GeV where the factor was the smallest. The pion-rejection factor for multi-hadron events is indicated in 
Figure by the dashed line. The average value for multi-hadron events was 8600. In the real analysis we 
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Figure 21: Various distributions for electron candidates: (a) the distance between a TPC track and its 
closest BCL cluster, (b) E/P, (c) the confidence level (CL) for the TOF, and (d) dE/dx. Each of these 
distributions was obtained with all of the cuts except for that on the plotted quantity. 
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Figure 22: dE/dx distributions for each Py-bin indicated in the figure. The smooth lines are the best-fit 
curves described in the text. 



29 




0.4 0.5 0.6 0,7 0.80.9 1 -_. 2 , ~ .3 T . 

electron P T (GeV) 

Figure 23: Detection efficiency of electron tracks for single-track events. 
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Figure 24: Pion-rejection factor and electron efficiency for two-photon hadronic events; the dashed curve 
is the pion-rejection factor (right scale) and the solid one is the electron efficiency (left scale). 
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counted the number of electrons by fitting the dE/dx distributions. Therefore, in fact, the pion-rejection 
factor was far better than these values. 

10 Summary 

We have developed a method to identify electrons inside hadronic jets using the time-projection chamber 
and the lead-glass calorimeter in the TOPAZ detector system at TRISTAN. The energy-loss measurement 
with the time-projection chamber and the energy measurement with the lead-glass calorimeter provided 
good electron identification over a wide momentum range. In particular, low- Pt {Pt > 0.4 GeV) electron 
identification was proven to be possible. This is powerful for studying the open-charm production in 
two-photon processes. The pion-rejection factors and the electron efficiencies inside hadronic jets were 
obtained to be 163 and 68.4% for high-P T electrons and 137 and 42.7% for low-P T electrons in singlc- 
photon-exchange process, and 8600 and 36.0% for two-photon processes. 
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